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Gap solitons in a model of a hollow optical fiber 
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We introduce a models for two coupled waves propagating in a hollow-core fiber: a linear dis- 
persionless core mode, and a dispersive nonlinear quasi-surface one. The linear coupling between 
them may open a bandgap, through the mechanism of the avoidance of crossing between dispersion 
curves. The third-order dispersion of the quasi-surface mode is necessary for the existence of the 
gap. Numerical investigation reveals that the entire bandgap is filled with solitons, and they all 
are stable in direct simulations. The gap-soliton (GS) family is extended to include pulses moving 
relative to the given reference frame, up to limit values of the corresponding boost S, beyond which 
the solitons do not exists. The limit values are nonsymmetric for S > and S < 0. The extended 
gap is also entirely filled with the GSs, all of which are stable in simulations. Recently observed 
solitons in hollow-core photonic-crystal fibers may belong to this GS family. 



The concept of gap solitons (GSs), that was introduced 
by Mills et al. (l| , has drawn a great deal of interest in 
optics and other fields where nonlinear waves are impor- 
tant. After a family of exact solutions for the GSs was 
found in the model of a nonlinear optical fiber carrying 
a Bragg grating (BG) 0], they had been created in the 
experiment [3j , see further details in reviews 3. Optical 
GSs were also predicted in various BG configurations in 
the spatial domain and, very recently, solitons of the 
gap type have been observed in a discrete array of optical 
waveguides and in Bose-Einstein condensates |3|- 

A new optical medium has recently become available 
in the form of photonic-crystal fibers (PCFs), with a 
crystal-like structure of voids running parallel to the 
fiber's axis |||. An important variety of the PCF is a 
hollow-core one, which provides a unique means for de- 
livery of extremely powerful light signals, with the energy 
of up to 100 pJ, over the distance of up to 200 m j^] - 
^ll|. The hollow-core PCF may also be a basic element 
of soliton-laser schemes the general case, theo- 

retical models for ultra-short solitons in a hollow PCF 
are complicated, see, e.g. , Ref. J^. However, in some 
cases the experiment [l4j and detailed numerical calcu- 
lations 10, 01 suggest to use an approximation which 
includes a linear dispersionless mode, v, propagating in 
the gas-filled core, and a quasi-surface one, u, that pene- 
trates into the silica and may therefore be dispersive and 
nonlinear. Another medium in which a similar interplay 
between the two modes may be expected, is a hollow- 
core multi-layer fiber; gap-like ("cutoff') solitons have 
been recently predicted in it [T^ . In this work, we aim 
to demonstrate that the coupled-mode system governing 
the co-propagation of the core and quasi-surface waves 
may open a spectral bandgap, which gives rise to a two- 
parameter family of stable GSs. The family is a novel 
one, being quite different from the well-studied standard 
example Q. 

A straightforward derivation, foUowingthe well-known 
principles of the coupled-mode theory [j, 0|, makes it 
possible to cast coupled equations for local amplitudes of 
the above-mentioned dispersive and dispersionless modes 



in the following normalized form (details of the derivation 
will be presented elsewhere): 

iUz — icUr + -;^UrT + il^T-TT + + V = 0, (1) 

iVz + icVr + u = 0. (2) 

Here, z and r are, as usual, the propagation distance and 
reduced time for the unidirectional propagation, the ef- 
fective coefficient of the Kerr nonlinearity in the u-mode 
and inter-mode linear coupling constant are scaled to be 
1, the second-order dispersion of the u-mode is assumed 
anomalous, and its coefficient is also normalized to be 
1, and 2c is the group- velocity mismatch between the 
two modes [the respective coefficients in Eqs. (QJ and 
are made symmetric (t*c) by adjusting the refer- 
ence frame; in fact, a group- velocity asymmetry will be 
re- introduced below, when looking for moving solitons]. 
Note that phase-velocity terms can be eliminated, and 
the coefficients in front of iuz and ivz can be made equal 
by means of obvious linear transformations, which is im- 
plied in the equations. 

Because the solitons observed in the hollow-core PCFs 
are very short, with the temporal width ~ 200 fs, and the 
carrier wavelength may be close to the zero-dispersion 
point IjI, the model must take into regard the third- 
order dispersion - ^JIj which is accounted for by the 
coefficient 7 in Eq. Thus, the normalized system 

contains two free parameters, c and 7. 

Without displaying details of the derivation, we note 
that Az = 1 and Ar = 1 in the normalized units 
correspond, typically, to 10 cm and 100 fs in the real- 
world PCFs, and the nonlinear coefficient in Eq. QJ, set 
equal to 1, may actually range between 10"^ and 10"^ 
W~^cm~^ in physical units |9[ll0l|. 

An effect that the model does not include is the soli- 
ton's self-frequency shift due to the intra-pulse stimu- 
lated Raman scattering (generally, it is an essential effect 
for ultra-short solitons flTll). From the experiment, it is 
known that this effect may be both significant and negli- 
gible in the hollow-core PCFs, depending, in particular, 
on the choice of the core-filling gas 0|. In this Letter, 
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we focus on the simplest conservative model which does 
not include the self- frequency shift; its role will be con- 
sidered in detail elsewhere. We do not include either the 
fourth-order dispersion, self-steepening, intrinsic nonlin- 
earity of the core mode, and loss. Estimates demonstrate 
that these factors are not expected to change the results 
dramatically; in particular, the net loss accumulated over 
the experimentally relevant transmission length is ~ 1 dB 

UM 

First, it is necessary to find the spectrum of the lin- 
earized version of the system. Substituting (m, v) ^ 
exp (ikx — icor), we obtain a dispersion relation that can 
be easily solved, yielding two branches, 

4fci,2 = -a;2(l-H27w)± (3) 
^/^^(l + 27w)2 + 8ccj2(2c + UJ + 2juj^) + 16. 

Further consideration demonstrates that, once we set 
c > by definition, the spectrum Q gives rise to a 
bandgap, where solitons may exist, only in the case of 
7 > 0, which we assume below (this case is quite possi- 
ble, as PCFs allow one to engineer the dispersion of the 
guided modes). The formation of the bandgap is actually 
stipulated by the principle of the avoidance of crossing 
between dispersion curves of the decoupled equations (Q) 
and (|2Jl; it is known that this principle is a crucial factor 
in the explanation of dispersion properties of the PCFs 

MM 

A typical example of the spectrum that contains the 
bandgap is shown in Fig. ^ The gap-existence region 
covers nearly the entire quadrant (c > 0, 7 > 0) in the 
system's parameter plane. We also note that the model 
never gives rise to more than one gap. An asymptotic ex- 
pansion of Eq. shows that, for small c, the minimum 
value of the third-order-dispersion coefficient, which is 
necessary for the existence of the bandgap, scales so that 

/mm ^ 

In the gap, stationary soliton solution were 
sought for by the substitution of t), r)} — 

{U{T]k),V{T]k)}e:xjp{ikz) in Eqs. and lO and 

solving the resulting ODEs by means of the relaxation 
numerical method. Note that the functions U and V in 
this ansatz cannot be assumed real; however, they obey 
symmetry restrictions, so that Rc {U{t)) and Re(V'(r)) 
are even functions of r, while Im([/(r)) and Im(V"(r)) 
are odd. Stability of the GS solutions was then tested in 
direct simulations. Additionally, it was tested by means 
of the Vakhitov-Kolokolov (VK) criterion ^T§\: for the 
soliton family, the condition for the stability against 
perturbations with real eigenvalues is dE/dk > 0, where 
the energy is 



m = 



-\-oc 



\U{T;k)\' + \V{T;k)f 



dr. 



(4) 



As a result, it has been concluded that the entire 
bandgap is filled by the solitons, and they all are stable 
in direct simulations (extensive scanning of the parame- 
ter space has not turned up any example of an unstable 



1 - dispersion curve k(M) 

2 - band gap 

3 - soliton energy as function of k: E(k) 




FIG. 1: The solid curves are two branches of the 
dispersion relation for c = 1, 7 = 0.3. The 

bandgap is located between the branches. The bold 
dashed curve inside the bandgap shows the dependence 
E{k) for the soliton family found in the gap. As is 
seen, the family fills the entire bandgap, and it al- 
ways satisfies the VK stability condition, dE/dk > 0. 



soliton). Besides that, all the GSs satisfy the VK crite- 
rion, see Fig. A typical shape of the stationary GS is 
displayed in Fig. 12 and its stability against finite per- 
turbations is illustrated by Fig. El 

A natural extension of the GS family is generated by 
lending the soliton a finite velocity, in the reference frame 
of Eqs. 1^ and Note that, as well as the ordinary 
GS model, the present one features no Galilean (or other) 
invariance that would generate moving solitons automat- 
ically. To develop this extension, we rewrite Eqs. 
and lO in the boosted reference frame, replacing the in- 
dependent variables (z, r) by {z,t' = t — S z), where 5 is 
the velocity shift. The transformation replaces the terms 
i {uz, Vz) by I {uz — 5 Ur,Vz — 5 Vt) and, accordingly, in 
the dispersion relation fc is replaced by fc -f J w. As a 
result, the bandgap changes its shape with the increase 
of 5. 

Naturally, the gap closes when \5\ is too large, and 
it is necessary to find the corresponding maximum and 
minimum values of the velocity shift, (5max and (5min, up 
to which the bandgap exists. Straightforward analysis 
of the dispersion relation demonstrates that Jmax — c, 
which has a simple explanation: the above transforma- 
tion replaces the coefficient c in Eq. ^ by c — (5, and the 
gap closes down when this combination vanishes. The 
other limit value can be found numerically from the cor- 
responding algebra. The final result is presented in Fig. 
01 that shows (5,„ax and 5mm as functions of c for given 
7. A noteworthy feature of this diagram is that the gap- 
supporting region continues, for 5 ^ 0, to negative c. 

Solving the stationary equations in the boosted refer- 
ence frame numerically, we have concluded that the entire 
gap is again completely filled with solitons. The boosted 
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FIG. 3: Stable evolution of the soliton from 
Fig. 12 to which an arbitrary initial perturba- 
tion was added, whose amplitude is 2% of the soli- 
ton's amplitude. The evolution is of the field 
|w(2:,t)| is similar to that displayed here for |u(z, r)|. 
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(b) 



FIG. 2: A generic example of the gap soliton (with k = 
and (5 = 0) found inside the bandgap shown in Fig. ^ 
Note that the absolute values of the fields, |C/(t)| and 
|y(T)|, always have an ordinary single- humped shape. 



-0.5 0.5 1 1.5 2 

soliton vetocity (c) 

FIG. 4: In the reference frame boosted by the velocity 
shift S, the bandgap exists in the interval < S < 
(5niax = c. The two lines show and (5max vs. c. 



solitons obey the same symmetry constraints as in the 
case of (5 = 0, so that the real and imaginary parts of the 
stationary fields are, respectively, even and odd functions 
of t' . The fields \u{z,t')\ and |w(z,t')| are always single- 
humped. A characteristic example of the boosted soliton 
in shown in Fig. 5. Finally, adding arbitrary initial per- 
turbations shows that all the moving solitons are stable, 
in direct simulations, everywhere in the bandgap. 

In conclusion, we have proposed the simplest model for 
two co-propagating waves in the hollow-core fiber, one a 
linear dispersionless core mode, and the other a disper- 
sive nonlinear quasi-surface mode. The linear coupling 
between them may open a bandgap, through the mech- 
anism of the crossing avoidance. The third-order disper- 



sion of the latter mode is a necessary condition for the 
existence of the gap. The numerical results demonstrate 
that the entire bandgap is filled with solitons, and they 
all are stable in direct simulations. The gap-soliton (GS) 
family has been extended to include the boosted solu- 
tions, up to the limit values of the boost S (these values 
are different for positive and negative S). 

Solitons which were recently observed in hollow-core 
photonic-crystal fibers may belong to this GS family. 
The theoretical results suggest to extend the experimen- 
tal study of these solitons, and of the "cutoff' solitons 
in hollow-core multilayer fibers, in order to identify their 
plausibly GS nature. On the other hand, the theoreti- 
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cal model needs to be extended by inclusion of the self- 
frequency shift and some other features. Results of the 
extension will be reported elsewhere. 
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FIG. 5: A generic example of the boosted gap soli- 
ton with S = —0.5, found at c = and 7 — 0.3. 
Note that, for c = 0, the solitons with (5 = do 
not exist, see Fig. ^ In terms of \u\ and \v\, the 
boosted solitons always have a single-humped shape. 



After the submission of this paper for publication, a 
work by D.V. Skryabin, Opt. Exp. 12, 4841 (2004) has 
appeared, which considers a very similar model and also 
reports finding stable solitons in it. 
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